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The solution of the classic problem of stress in a rotating elastic disk or cylinder, as solved in 
, standard texts on elasticity theory, has two features: dynamical equations are used that are valid 

only in an inertial frame of reference, and quadratic terms are dropped in displacement gradient 
in the definition of the strain. I show that, in an inertial frame of reference where the dynamical 
equations are valid, it is incorrect to drop the quadratic terms because they are as large as the 
linear terms that are kept. I provide an alternate formulation of the problem by transforming the 
dynamical equations to a corotating frame of reference of the disk/cylinder, where dropping the 
quadratic terms in displacement gradient is justified. The analysis shows that the classic textbook 
. derivation of stress and strain must be interpreted as being carried out in the corotating frame of 

the medium. 
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I. INTRODUCTION 

The problem of stresses in rotating disks and cylinders is important in practical applications to rotating machinery, 
such as turbines and generators, and wherever large rotational speeds are used. The textbook problem of stresses 
■ in elastic rotating disks and cylinders, using the assumption of plane strain or plane stress, is published in classic 
texts, such as Love jj| , Landau and Lifshitz g| , Nadai || , Sechler Q , Timoshenko and Goodier J5| , and Volterra and 
Gaines [§. The standard approach presented in these texts has two characteristic features: 

i-C ' 1. Newton's second law of motion is applied in an inertial frame of reference to derive dynamical equations for the 

continuum (see Eq. (TP) below), and 
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2. quadratic terms in displacement gradient are dropped in the definition of the strain tensor (see Eq.(EOj) below) 



In this paper, I show that, for a rotating elastic body, the second feature of the solution is inconsistent with the first: 
dropping the quadratic terms in the displacement gradient is an unjustified approximation in an inertial frame of 
reference. In what follows, I refer to the method that is employed in Ref. [|IJ, @, 0, EL pL 0] as 'the standard method', 
t-H . and for brevity, I will refer to a cylinder as a generalization of both a disk and a cylinder. 

The classic problem of stress in an elastic rotating cylinder is complex because the undeformed reference state 
of the body is the non-rotating state. The deformed state is one of steady-state rotation. The analysis of the 
problem must connect the non-rotating reference state to the rotating stressed/strained state. These two states are 



^ , typically connected by large angles of rotation. When large angles of rotation are present, the quadratic terms in the 
. ^ ■ displacement gradient cannot be dropped (in an inertial frame of reference) in the definition of the strain Jt], || ^ . The 
problem of stress analysis when large-angle rotations jare present is well known and has been discussed by a number 



of authors in general contexts, see for example ||, However, large-angle rotations in the problem of a rotating 
elastic cylinder have not been dealt with in a technically correct manner, because quadratic strain gradient terms are 
incorrectly dropped in the 'standard method' [Q, || . 

In this work, I formulate the elastic problem of a rotating cylinder in a frame of reference that is corotating with 
the material. In this corotating frame, the quadratic terms in the displacement gradient can be dropped, and the 
5^ ■ resulting differential equations are linear and can be solved. 

In section II, I review the 'standard method' of solution used in Ref. |l], |[ [| f§, ||j H an d show that for a rotating 
cylinder the displacement gradient in an inertial frame of reference is of order unity, and therefore quadratic terms 
(in strain tensor definition) cannot be dropped when compared to the linear terms. Section III contains the bulk 
of the analysis. I describe the corotating systems of coordinates and the transformation of the velocity field to the 
corotating frame. I use the velocity transformation rules to transform the dynamical Eq. ([j]) from the inertial frame 
to the corotating frame (see Eq. ( J6l| ) or (|6^)), where extra terms arise known as the centrifugal acceleration and the 
coriolis acceleration. In section IV, I write the explicit component equations for stress (in cylindrical coordinates) 
for the rotating elastic cylinder in its corotating frame. To display the resulting solution concretely, I derive the 
well-known formula for the stress in the rotating cylinder for the case of plane stress, as computed in the corotating 
frame. Stress is an objective tensor, i.e., stress is independent of observer motion jl0| [n]], so the physical meaning 
of stress in the corotating frame is the same as in the inertial frame. Therefore, the stress field components in the 
corotating frame are equal to the stress field components in the inertial frame, see Eq. (p6|). 
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II. STANDARD SOLUTION METHOD 



In the 'standard method' |jL ^, ||, 0,0,0], the stress analysis of elastic rotating cylinders starts with the dynamical 
equations, which, in generalized curvilinear coordinates are given by [|[ |o[ [ll| 

a%+pf k =pa k (1) 

where a k ^ are the contravariant components of the stress tensor, f k is the vector body force, and a k is the acceleration 
vector. In Eq. dll), repeated indices are summed and the semicolon indicates covariant differentiation with respect to 



the coordinates. Expressed in terms of the velocity field in spatial coordinates, the acceleration is given by fid, 11 



«' ^ ""^ (2) 

where v a is the velocity field, and the semicolon indicates covariant differentiation with respect to the coordinates, 
and v^v k .^ is called the convectivc term. In Eq. (Q), the stress a k ^ , acceleration a k , and body force f k , are generally 
time dependent. Equation (|l]) is derived by applying Newton's second law of motion to an element of the medium. 
Newton's second law is valid only in an inertial frame of reference, and consequently the validity of Eq. ([!]) is limited 
to inertial frames of reference. 

In the 'standard method' of solution, Eq. (jl]) is applied by invoking an "effective body force" , of magnitude equal to 
the centrifugal force in the rotating frame. In the inertial frame, there is actually no effective force (such as Coriolis 
or cetrifugal force). For the case of a body rotating about its principle axis, a more careful determination of the terms 
f k — a k in Eq. (|]) comes from setting the body force to zero (or setting equal to some applied force) and computing 
the material acceleration a k for a given body motion. For a rigid body, or a uniform density elastic cylinder that is 
rotating about its axis of symmetry at a constant angular velocity w , the Cartesian velocity field components are: 
v = —uj y, v 2 = uj x, and v 3 — 0, where superscripts 1,2,3 indicate components on the Cartesian basis vectors 
associated with the x,y,z-axes (in the inertial frame). Corresponding to this velocity field, the cylindrical components 
of the acceleration field are given by 

a* = ^+« b « fe ;6 =(-r^0,0) (3) 

where I have chosen the z-axis as the symmetry axis and the bar over the components indicates that they are in the 
inertial frame of reference in cylindrical coordinates. For the case where there are no body forces, with the acceleration 
in Eq. (||), Eq. (|j) in cylindrical coordinates leads to the three equations 

li 



13 i a -22 „„, ,2 



3 i 



ra 2A = -pro; 2 (4) 



a 12 1+ a 22 2 + a 23 3 + h 12 = (5) 

^13 l+(f 23 2 + (f 33 3 + ^! = (6) 

where the superscripts 1,2,3 enumerate tensor components on the r,4>,z coordinate basis vectors respectively, in 
cylindrical coordinates and the commas indicate partial differentiation with respect to these coordinates. 

For steady rotation at a uniform angular velocity oj 0i and assuming the absence of elastic waves, there is rotational 
symmetry about the z-axis so the stress components do not depend on azimuthal angle (f>. Therefore, all derivatives 
with respect to <j) are zero, leading to the equations: 

- n .i + - 13 3 + — - ^ 22 = -P^l (7) 
r 

-12 , -23 , ^ -12 n /n\ 

r 

^ 13 i + ^ 33 3 + ^ - (9) 

rrj> 



I introduce physical components of stress, a rr , cr^cr 22 , <7 r< * ) , cr rz , and a^ z , with units of force per unit area and 
which are related to the tensor components a 11 , a 22 , ct 33 , a 12 , <r 13 , and ct 23 , by |ll| 

a rr = a 11 , (10) 
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d** = r 2 a 22 (11) 

a zz = a 33 (12) 

a r * = ra 12 (13) 

a rz = a 13 (14) 

= ra 23 (15) 

Expressing Eq. (ff|)-(B) in terms of the physical components, I obtain the well-known equations valid in an inertial 
frame of reference ^fjfl, ||, [|, ||, |(| , 

da rr da rz a rr - , 

— + — + : = ~Pr«t (16) 



z r 



= (18) 



9r \ r / r dz r 2 



9r <9z r 

Note that Eqs. (Jl6|)-(|T8"|) have been derived using Newtons's second law, and so they are valid only in an inertial 
frame of reference. In particular, Eqs. (|l6|)-(|l8|) are not valid in a rotating frame of reference. 

When a rotating disk or cylinder is analyzed, the assumption of plane stress or plane strain is often made. In 
both cases, stresses must be related to strains by constitutive equations. For the simplest case of a homogeneous, 
isotropic, perfectly elastic body, the constitutive equations in curvilinear coordinates in an inertial frame can be 
written as g [HJ |ll| 

a tk = \eg ik + 2 M e ifc (19) 

where A and \i are the Lame material constants, e lk are the contravariant strain tensor components, e = e a a is the 
contraction of the strain tensor, and g zk are the contravariant metric tensor components. 
The Eulerian strain tensor e,fc is related to the displacement field u 1 by |l(| |ll|] 

e-ik = ^ { u i-,k + u k , % + u m;i u m . k ) (20) 

In the 'standard method' of solution |], |, §, |, §, the quadratic terms iifji-iU Eire dropped, which leads to linear 
equations that can be solved (for example, by using the Airy stress function Jl2]|). 

However, dropping the quadratic terms in Eq. (g0|) is not justified for a rotating body because these (dimensionless) 
terms u m -,i are of order unity. To prove this assertion, it is sufficient to consider the limiting case of a rigid body 
in steady-state rotation at constant angular speed u> - The deformation mapping function gives the coordinates z k 
(here taken to be Cartesian) of a particle at time t in terms of the particle's coordinates Z k in some reference state 
(configuration) at time t — t a : 

z k = z k (Z m ,t) (21) 

so that z k (Z m ,i Q ) = Z k . The deformation mapping function has an inverse, which I quote here for later reference 

Z m = Z m (z k ,t) (22) 

Both z k and Z m refer to the same Cartesian coordinate system. The coordinates of a particle initially at Z k at 
t = t Q = rotating about the z-axis are given by the deformation mapping function 



R k m Z m (23) 



where the orthogonal matrix R m is given by 



cos ui t sin oj t 

R k m = \ -smu> t cosuj t | (24) 
1 



The displacement vector field for this deformation mapping function is given by 1 13 

u = u m I m = (z m - Z m )I m = (51" - R m k )z k l m (25) 
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where I m are the unit Cartesian basis vectors and R m k is the transpose matrix that satisfies 

R m k R k l =Sr (26) 

where Sf l = +1 if m = / and if m ^ I. Therefore, from Eq. (|25|), it is clear that gradients of displacement u m -k 
appearing in Eq ( ^o|) are of order unity and therefore the quadratic terms u rn ^u m . k cannot be dropped because they are 
not small. More specifically, the dropped terms in Eq. (|20| ) vary in time between -2 and (in Cartesian components): 

1 / -l + cosui t 0\ 

-u m;l u m . fc = -1 + cos^t (27) 

1 \ 0/ 

The above calculation was done for a rigid body, but clearly, a similar error is introduced for elastic bodies. Therefore, 
in general, for a rotating elastic body, the quadratic terms in displacement gradients in Eq. (|2p|) cannot be dropped . 

In cylindrical components, the relation between the physical components (see Ref. [pLOj, of strain, e rr and e^, 
and physical components of the displacement field, (u r ,u^,U z ), is given by 

du r 1 

Crr = 57 ~ 2 

eu = [ (U " )2 + (29) 

(I use a bar over e rr and to indicate that these are cylindrical components, and indices rr and <jxj> (as distict from 
11 and 22) to indicate that these are physical components and not tensor components. See the Appendix and Table 
I and II for notation conventions.) In Eq. ( p8| ) and (p9|), I have assumed that there is no dependence on </> and z, so I 
have set derivatives with respect to these variables to zero. 

In the 'standard method' of solving for the stress in a rotating cylinder [Q 0, ||, 0, |L El, the quadratic terms in 
Eq. ( p8| ) and (|2^) are incorrectly dropped. 

The straight forward approach to correctly studying the stresses in a rotating disk or cylinder, involves keeping the 
quadratic terms in displacement gradient in Eq. (|2(]). However, this approach does not appear promising because it 
leads to insoluble nonlinear differential equations. In the next section, I approach the problem by using a transforma- 
tion to a corotating frame of reference, in which dropping the quadratic terms can be justified for moderate angular 
velocity of rotation lo . 



du r \ ( du^\ ( du z 
dr J \ dr J \ dr 



(28) 



III. TRANSFORMATION TO THE ROTATING FRAME 



As discussed in the introduction, the problem of an elastic rotating cylinder is complicated because the un- 
stressed/unstrained reference state is the non-rotating state, while the stressed (strained) state is rotating, and these 
two states are typically related by a large (time-dependent) angle. The analysis of the rotating disk or cylinder must 
relate the stresses in the rotating state to the reference configuration, which I take to be the non-rotating state. I 
define a transformation from an inertial frame of reference to the corotating frame of reference of the cylinder. This 
transformation provides a relation between the rotating stressed state and the non-rotating reference configuration. 



A. Coordinate Systems 



Starting from an inertial frame of reference, S, defined by the Cartesian coordinates z k = (x,y,z), I make a 
transformation to a rotating frame of reference S' . The rotating frame will be corotating with the cylinder so that 
in this frame S' the azimuthal velocity field will be zero at all times. The transformation from an inertial system of 
coordinates to a rotating system of coordinates is most simply done using Cartesian coordinates. On the other hand, 
the assumed cylindrical symmetry of the problem begs for use of cylindrical coordinates. Hence I will make use of four 
systems of coordinates. In the inertial frame S, I have two systems of coordinates: a Cartesian system of coordinates 
z k = Zk = (x,y,z), and a cylindrical system of coordinates x l — (r, <j>, z). In the corotating frame of reference, 
S", I have a Cartesian system of coordinates z' k — (x',y',z') and a cylindrical coordinate system x n = (r' , <j>' , z'). 
These coordinates are summarized in Table | and the Appendix. I also introduce notation for tensor components 
in each of the four coordinate systems, see Table [II] . The Cartesian components of the stress tensor in the inertial 
frame S will be denoted by a %k . In the same inertial frame S, the cylindrical components of stress will be a lk . The 
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Cartesian components of stress in the corotating frame S' will have a prime, a nk . In this same corotating frame, S', 
the cylindrical components of stress will be denoted by using a tilde, a lk . 

From the vantage point of an inertial frame of reference, S, with Cartesian coordinates z k , consider a cylinder 
whose symmetry axis is aligned and colocated with the coordinate z-axis. At time t = — oo, take the cylinder to 
be non-rotating. Now assume that in the distant past, around the time t ~ — T, the cylinder begins a slow angular 
acceleration lasting a long time, on the order of 1/e, where 1/e << T. An example of such an angular acceleration 
function is 

w(t) = ^u [1 + tanh(e(> + T))] (30) 

where I assume that r << 1/e << T and r is the longest time constant in the problem. This inequality states that 
the acceleration occurs slowly, r << 1/e, slower than any time scale in the problem, and that this acceleration occurs 
in the distant past, 1/e << T, so that at t = 0, I have a steady-state situation of a cylinder rotating at constant 
angular speed oj . By slowly accelerating the cylinder, I avoid introducing modes of vibration. As the cylinder's 
angular velocity increases from t — — oo, each particle comprising the cylinder moves along a spiral trajectory (with 
increasing radius). From the point of view of the inertial frame S, the stresses on a given element of the medium 
(particle) are such that they cause the particle to experience an acceleration, moving along the spiral path. At t = 0, 
the cylinder has achieved its maximum angular velocity u) . Due to the assumption of a perfectly elastic medium, at 
t = the velocity field has zero radial component; all particles of the cylinder are moving azimuthally (in a plane 
perpendicular to the z-axis with zero radial component). The velocity field is that of a rigid body and the acceleration 
field is given by Eq. (|3|). 

Now I introduce the corotating frame of reference, S', with the Cartesian coordinates z' k , whose angular velocity of 
rotation is equal to that of the cylinder at all times. The coordinates z k (= z k ) and z' k (= z ) are related by 

z'^A lk (t)z k (31) 

where the time dependent matrix Ai k (t) is given by 

cos(0 - O ) sin(0 - 

A tk (t)= ( -sin(0-0 o ) coa(e-9 ) | (32) 




where is a function of time given by the integral of tu(t): 




= ^ ^ + T ) + lo § cosh(et + e T) + log 2] (33) 

and 0(0) = O . 

By construction, in the corotating frame S' the particles comprising the material are not rotating about the z'-axis; 
there is zero azimuthal component of the velocity field at all times. As the angular velocity u>(t) increases from 
t = — oo, each particle comprising the cylinder experiences an increasing effective centrifugal force that displaces the 
particle to a larger radius. In this rotating frame S', there will (in general) also be a Coriolis force. However, in S', 
for moderate angular speed u> , the strain will be small, and the gradients of the displacement field will also be 
small. Consequently, dropping the quadratic terms u m -iU rn . k in Eq. ( p0| ) will provide a good approximation to e^fc- 

Note that the transformation that relates cylindrical components in inertial frame S and rotating frame S' is given 
by the identity matrix 

dx l dx l dz a dz' m ■ 

dx 71 = dz^ dz' m dx' k k ( ' 

Furthermore, the time dependent transformation between inertial cylindrical coordinates x % — (r, <f>, z) in S and 
corotating cylindrical coordinates x' % = (r', 0', z') in 5", is given by 

r = r 

<f>' = 0-(0(t)-0 o ) (35) 

z' = z 

where 0(t) is given by Eq. (|33|). The relation between cylindrical stress components <j %k in the inertial frame S and 
cylindrical stress components a lk in the corotating frame S' is 

* ik ^ = §fa^* ab (x' n )=* ik (x' n ) (36) 
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Of course the transformation in Eq. ( |36| ) must be used so that the components are referring to the same physical 
point in space having coordinates x n = (r, 0, z) and x' n — (r', $ , z'), where the time-dependent relation between x n 
and x' n is given by Eq. fl35|). 

The Eulerian strain tensor transforms in a more complicated manner. A general deformation is given by 

x l =x l {X k ,t) (37) 

where a particle at time t = t a in the reference configuration has (curvilinear) coordinates X k , and in the deformed 
state at time t the particle has coordinates x l (in the same curvilinear coordinate system). The Eulerian strain tensor 
eij(X, x) depends on two points: X in the reference configuration and x in the deformed state. Consequently, under a 
general coordinate transformation to a moving frame, x % — > x l% = h % (x k , t), the Eulerian strain is a two-point tensor, 
which transforms as a second rank tensor under transformation of deformed coordinates 

x l -> x 11 = h\x k ,t) (38) 

and transforms as a scalar under transformation of reference state coordinates 

X 1 -> X' 1 = h i (X k ,t a ) (39) 



so that H 



e mn (X,x) = e ik (X',x') = e mn (X',x') (40) 



where x and x', and X and X', are related by Eq. ( |38| ) and (|39|), and I used Eq. (|34|). Therefore, the cylindrical 
components of strain in the inertial frame S 7 e mn (X,x), are equal to the cylindrical components of strain in the 
rotating frame S', e.i k {X' ,x'). Finally, since I assume rotational symmetry about the z-axis (and z'-axis) so that all 
physical quantities have no dependence on </> or <!>, which are the azimuthal coordinates of the point in the deformed 
state, x = x l — (r, (j>, z), and coordinates of the point in the reference configuration, X = X 1 = (R, $, Z). Because of 
the nature of the transformation to the rotating frame in Eq. (|35|), Eq. ( p0| ) can be used with r' — r, R' — R, z' = z, 
and Z' = Z. Therefore, the tensor components in the corotating frame are identical to the components in the inertial 
frame. 



B. Transformation of Lagrangean and Eulerian Velocities 



The motion of a particle in the cylinder is given by the deformation mapping function. In the inertial frame S, using 
Cartesian coordinates, the motion of the particle is given by Eq. (|2l|), where the particle in the reference configuration 
at time t = t a has coordinates Z m . The coordinates Z m label the particle in the Lagrangean desription. Using the 
transformation to the rotating frame in Eq. (|3l|), the motion of the particle with label Z m with respect to the rotating 
S' frame Cartesian coordinates is given by 

z lk {Z m ,t)=A k] {t)zJ{Z m ,t) (41) 

In discussing the transformation to the rotating system of coordinates, I must distinguish between the velocity of a 
given particle in the medium (the Lagrangean picture) and the velocity field (the Eulerian picture). The Lagrangean 
velocity Vi(S; Z m ;t) of a given particle (whose coordinates in the reference configuration are Z m ) with respect to the 
inertial frame S is defined as the partial time derivative of that particle's z 4 -coordinates, when holding Z m constant 

v,(5;^;*) = ^^ (42) 

The Eulerian velocity field, v.i(z k ,t), with respect to the frame S is a function of coordinates z k and time t and is 
related to the Lagrangean (particle) velocity by 

v t (S;Z m ;t) = v i (S;Z m (z k ;t),t) = Vl (z k ,t) (43) 

where I used Eq. (|2^) to express the particle coordinate Z m in terms of its position z k at time t. 

I can describe the same particle's velocity (whose coordinates in the (inertial frame) reference configuration are 
Z m ) with respect to the rotating frame of reference S'. The velocity of this particle with respect to the rotating frame 
S' is 

v' i {S';Z m ;t) = A „ ' j = - [A ik (t)z k (Z m ,t)\ = A ik (t) z k (Z m ,t) + A ik (t)v k (S; Z m ;t) (44) 
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where I used Eq. ( |3l| ) to express the particle's coordinates in terms of coordinates in the S frame and the dot on A, Lk 
indicates differentiation with respect to time. 

The components Vi(S; Z m ;t) and v^(S"; Z m ;t) represent physically distinct vectors (geometric objects). Each of 
these vectors can be expressed on the other basis. In particular, according to the standard transformation rules for 
vector components, I have 

v' i (S';Z m ;t) = A lm v m (S'-Z m ;t) (45) 
v' m (S;Z m ;t) = A mi Yi(S;Z m ;t) (46) 



Equations ([15|) and (46) are the standard tensor transformation rules for vector components under the coordinate 
transformation given in Eq. (J3l|). In summary, I must distinguish between four (Cartesian component) velocities [18j: 



Vj(5; Z m ; t) = ii(Z m , t) = components on z-axes of particle velocity with respect to S 

Vi(S'; Z m ;t) = components on z-axes of particle velocity with respect to S' 

Vj-(iS"; Z m ;t) — z[{Z m , t) — components on z'-axes of particle velocity with respect to 5" 

v' m (S;Z m ;t) = components on z'-axes of particle velocity with respect to S 

These four velocities are related. Multiplying Eq. (^4|) by A im (t), summing over index i and using the orthogonal 
matrix properties 

A in (t)A im (t) = 5 nm (47) 
A m (t)A mi {t) = 5 nm (48) 

leads to || 

v n {S'; Z m ; t) = v n (S; Z m ;t)+ u nk {S> \ S) z k (Z m ,t) (49) 

where uj^^S 1 , S) is the Cartesian angular velocity tensor of frame 5" with respect to frame S: 

u nk (S',S)=A in {t)A ik (t) (50) 

The tensor uj n k(S',S) describes the time-dependent rotation of frame S' with respect to frame S. Similarly, using 
Eq. (^l|), substituting the inverse relation 

Zi (Z m ,t)=A ni (t)z' n (Z m ,t) (51) 

carrying out the time differentiation, multiplying by A m i, summing over i and use of the orthogonality relations in 
Eq. ( jiq ) leads to 

v' n (S; Z m - 1) = v' n (S'; Z m ; t) + co' kn (S\ S) z' k (Z m ,t) (52) 
where the angular velocity tensor components are expressed with respect to the S' frame Cartesian basis: 

u>' nm (S',S) = AmiAni = A ni A mk oJik(S',S) (53) 



Equation (|49j) and (J52J) are the well-known rules for transformating particle velocity to a rotating frame of reference p8| . 

Next, I derive the equation that relates the Cartesian components of the velocity field in S, Vi(z k , t), to the velocity 
field in S', u^(z' fe ,t). Equation (43) relates the velocity field in the S frame to the Lagrangean (particle) velocity. 
Similarly, the velocity field with respect to the S' frame is given by 

v' i (z' k ,t)=v' i (S';Z m ;t) = = v^S'; Z m (z k , t); t) = y'^S'; Z m (A nk z' n , t); t) (54) 



where I used the inverse relation z k — A nk z' n . Using Eq. 



and (M) in the left and right most terms in Eq. (p4), 



I obtain a relation between the velocity fields in frames S and S' 

v' l (z' k ,t)=A tk (t)z k + A tk (t)v k {z n ,t) (55) 
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In Eq. (p5|), the coordinates z' k and z n are related by Eq. (fil]). Multiplying Eq. ( |55| ) by Ai m , summing over index, 
z, using the inverse transformation in Eq. ( |5l| ) and 

oJ jk (S', S) = A mj A nk u>' mn (S', S) (56) 

leads to 

Vj(z m ,t) = A lJ (t)v' l (z' nl t)-A kj (t)u;' kn (S' 1 S)z , n (57) 

Eq. ( |57j ) is the desired rule for transformation of the Eulerian velocity field from the inertial frame S to the rotating 
frame 5". Note that the right side of Eq. (|57]) depends only on S' frame coordinates z' n and the left side depends 
on S frame coordinates z m . Furthermore, the velocity field components on the left side of Eq. ( [57| ) are taken on 
the (Cartesian) inertial frame S basis vectors, and on the right side all components are expressed on the (Cartesian) 
rotating frame S' basis vectors. 



C. Dynamical Equation in the Rotating Frame 

In what follows, I transform the momentum balance Eq. (Q) to the corotating system of coordinates S". For 
simplicity, I do this transformation using Cartesian coordinates for both the inertial frame S and corotating frame 5". 
I use the transformation of the velocity field given in Eq. (|5^) to compute the terms that appear in Eq. (Q). Taking 
the gradient of the velocity in Eq. (|57j) 



— — A nk Aij — - — -p- A mk A n j u)' nm (S , S) (58) 

o z k z n 

where I used the chain rule for differentiation g|- — A mk -^r- since z k and z' m are related by Eq. (|4l|). Next, I 
compute the time derivative of the velocity that occurs in Eq. (Q) and I express the right side in terms of S' frame 
components and coordinates: 

dUi^ZjYii ^) A I \ I I A ^^TTl A It A • I I A I II I r r\\ 

Q-j_ = A m i v m + An LO nm z rn — — + A m i —gj A r i uj rn z n — Ami ^mn z n ~ An 0J nm LO ln Z m (5J) 

where I have omitted the frame labels of the angular velocity and the coordinate arguments in the velocity. The 
divergence of the stress transforms as a vector under the transformation to the rotating frame 

^ = A ai (60) 
a z k dz' b 

Substituting Eqs. (|57|)~(|60|) into the inertial-frame momentum balance Eq. ([!]) and simplifying, leads to the dynam- 
ical equation for the (Cartesian) stress tensor in the rotating frame S' 

lda[ k _ dv[ , dv[ { t ! duj' in \ , , ! 

p~dzl Vn dzl + \ u *™ w ™ ~ ~dT ) Zn ~ im Vm [ ' 

where a' ik are the stress components in Cartesian coordinates z' k , v' n is the Eulerian velocity field that depends on 
z' k and t, and the (Cartesian) components of the angular velocity tensor in S' are given by Eq. (^3|). In Eq. (|6l]), 
all repeated subscripts are summed. Note that all velocities that appear in Eq. (ph) refer to the corotating frame 
S' and that all tensor components are taken on the S' frame Cartesian basis vectors. The first two terms in Eq. 
( |6l| ) are the acceleration (including the convective term) as seen in the corotating system of coordinates. The third 
term uj' im u' mn z' n is the centrifugal acceleration. The fourth term, io' in z' n is the angular acceleration. The last term, 
—2 ui' in v' n is the Coriolis acceleration. 

The dynamical Eq. (|6|) is a tensor equation; the quantities a' kn , u>' „, and v' n , are Cartesian tensors. Under 
orthogonal transformations from one Cartesian system to another, Eq. (pi]) is covariant: it has the same form. The 
group of symmetry operations may be extended to transformations between curvilinear coordinates by writing Eq. 
( |6l| ) in a manifestly covariant form as: 

1 r)?? f Bui * \ ~ 

- a* fe ;fc = — + v n v\^ n + 2 u, n ' v n + (u m * u n m + C (62) 
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where the tilde over each tensor indicates that the components are taken on generalized curvilinear coordinate (such 
as cylindrical) basis vectors in the rotating frame of reference S' . Equation ( p2| ) is general; it is valid for all motions 
and all materials. Equation (^32|) expresses Newton's law for a continuous medium in an arbitrary rotating frame of 
reference that is defined by a general (Cartesian) angular velocity tensor u>ik(S', S) = Au(t) Aik(t), which relates the 
inertial frame S and rotating frame S', with coordinates related by Eq. (|3l]). The quantities £™ are the contravariant 
components of the position vector in curvilinear coordinates x' m , which are related to the Cartesian position vector 
components z lk by 

r = %u ** (63) 

The rotating frame curvilinear components of stress, velocity and angular velocity, a ab , v a , and u) a b , are related to 

Ty, v' k , and u' ip 



their rotating frame Cartesian components, e^ -, v' k , and u' iit by: 



b dx"> dx' b 

° = (64) 

flT ,a 

" a = v ' k (65) 

- b dz» dx' b , 

w o = T\ — — Ft — T DO 

a dx' a dz'i 13 K ' 

where z' k and x' a are the Cartesian and curvilinear coordinates in the rotating frame 5". Note that the partial 
derivative with respect to coordinates z' k in Eq. (Jfjl) has been replaced by a covariant derivative with respect to 
the curvilinear coordinates x' k in Eq . (|62| ). See Tables I and II for summary of the notation. Identification of the 
meaning of the various terms in Eq. (|62|), such as the Coriolis acceleration and centrifugal acceleration is clear from 
comparison with Eq. (|6l|). 



IV. ROTATING CYLINDER EQUATIONS IN COROTATING COORDINATES 



In order to obtain the explicit equations for a rotating cylinder from Eq. (|62|), I need to compute v n , uJ m n , and £™. 
As described previously, I take the reference configuration of the cylinder to be the stationary cylinder at t — — oo in 
the inertial frame S. I assume that the cylinder experiences a slow angular acceleration (such as given by Eq. (pQ)) 
that lasts approximately time At = — 1/e, and has peak magnitude at t = —T, with 1/e << T. Equation (|62[) is 
general; it is valid for all motions and all materials. In what follows, I restrict my remarks to a perfectly elastic 
cylinder. 

In the inertial frame S, as the cylinder increases its angular velocity, the material particles of the cylinder move 
outward in a spiral path, with some motion in the z-direction. In the inertial frame S, the contravariant components 
of the velocity field in cylindrical coordinates are 

« i = (« 1 ,w(t),t? 3 ) (67) 

where v 1 and v 3 are the radial and z-components of velocity and where the azimuthal component u>(t) is given 
by Eq. (p0|). I obtain the velocity components in cylindrical coordinates in the corotating frame S' as follows: 
first, transform v % to inertial frame Cartesian components v % using the standard vector transformation rule between 
Cartesian and cylindrical coordinates. Next, use Eq. ( p7j ) to transform the inertial frame Cartesian componets v % 
to the rotating S' frame Cartesian velocity field v' 1 . Finally, use the standard tensor transformation rules, between 
Cartesian and cylindrical coordinates (both in the rotating frame), to transform the velocity field from Cartesian 
components v n to cylindrical (rotating frame) components v n , where I made use of the angular velocity components 

/ -1 0\ 

uj ik (S', S) = A u A lk = ui(t) +1 = u4(£', S) = A im A kn w mn (S', S) (68) 

\ J 

Following this procedure, I obtain the velocity components in cylindrical coordinates in the corotating frame S' 

& = (v 1 , 0, v 3 ) (69) 



where the azimuthal component of velocity is zero for all time, by construction of the corotating frame S', as expected. 
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In the corotating frame S', there is particle motion around the time t rj — T. However, at t — the particles have 
reached their new (deformed) steady-state positions and motion has ceased; the velocity field is given by 



D fe = (0,0,0) (70) 

which is the velocity field of a rigid body. Since the velocity field is zero at t = 0, the Coriolis acceleration term in 
Eq. ( |62| ) does not contribute in steady-state rotation. 

Using Eq. (|66|) , the angular velocity tensor in cylindrical coordinates in the corotating frame S' is 



Lo(t)/r' 

-r'u{t) ) (71) 


At t = 0, the time derivative of the angular velocity tensor is zero. Using this fact, and Eq. ( |7C| ) and (fn|), the stress 
Eqs. ( |62| ) in the corotating frame are given by 

a" 1+ a 13 3 + ^-r'a 22 = -pr'^ (72) 
^ 12 i+^ 23 3 + ^ 12 = (73) 
^ 13 i+^ 33 3 + ^ = (74) 

Equations ([7^)-(|74]) are the equations satisfied by the stress tensor in cylindrical com po nents at t = 0, in steady-state 
rotation in the corotating frame S' . In terms of physical components (see Eq. (|l0|)-(^5[)), Eq. (|72|)-(|7^) become: 

= ~pM (75) 



da rr 


da rz 


dr' 


dz' 


~r<j> \ _ 


1 dd* z 


J ' r' dz' 



r 



d °" di " ,77) 



dr' 



Note that Eq. (|75|)-([77|) in the corotating frame S' have the same form as Eq. (|l^)-(|l8|) in the inertial frame S 



However, the key point is that the corotating frame Eq. (|75[)-(77) have a distinct advantage: the quadratic strain 
gradient terms in the definition of the strain in Eq. (^0|) can be dropped because they are small in the corotating 
frame S' . The same is not true in the inertial frame S. 

To proceed with the solution in the corotating frame, the constitutive Eq. ( p~Qj ) (in inertial frame S) must be 
transformed to the corotating frame S' by taking cylindrical components in the S frame and using Eq. ([54]) to obtain 
an expression of the same form as in Eq. ( |l^ ) but in the corotating frame S' . In this transformation, the Lame 
constants are treated as invariants. So the constitutive relations in the corotating frame S' are the same as in the 
inertial frame S. 

In the corotating frame S', dropping the quadratic terms in displacement gradients, which relate the displacement 
field to strain, is justified since in the rotating frame S' these terms can be considered small for moderate angular 
velocity. Therefore, the solution in the corotating frame S' can proceed in an analogous way to that of the 'standard 
method', but I have not made the (incorrect) approximation of dropping quadratic displacement gradient terms in 
the inertial frame. 

The stress tensor is objective, so the stress in the rotating frame has the same meaning as in the inertial frame, see 
Eq. (|3^). The boundary conditions on the stress tensor components in the corotating frame are the same as in the 
inertial frame, due to the objectivity of stress tensor. Alternatively, one can verify that the boundary conditions on 
the stress in the rotating frame are the same as in the inertial frame EQ] . 



V. PLANE STRESS SOLUTION 



The solution of the problem of stress in a rotating cylinder in the corotating frame of reference now follows. The 
solution in the corotating frame parallels the solution in the 'standard method' |l|, g, |3| Q, || |fj, except that the 
incorrect approximation of dropping the quadratic strain gradient terms in the inertial frame is avoided. 
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I assume that in the inertial frame S, the cylinder is rotating at angular velocity uj and has radius b. Under the 
assumption of plane stress M , where 



~zz ~d>z ~rz n 

a = <j = <j = U 
with boundary condition of zero stress on the long peripheral surface: 



(78) 



(79) 



As mentioned above, stress is an objective tensor |TTj, so the physical meaning of the boundary conditions in 
Eq. ( fZ8| ) and ( |79] ) in the corotating frame S' are the same as the physical meaning of the analogous conditions in the 
inertial frame (as used, for example by Timoshenko ||). Also, as shown above, I may take the (transformed) linear 
elastic relations in the corotating frame to be of the same form as in the inertial frame: 



e rr = — (a rr -vo* 



E 

where Young's modulus E and Poisson's ratio v are related to the Lame constants by 



E 



(80) 
(81) 

(82) 
(83) 



^(3A + 2^) 
X + fi 
A 

2 (A + fj.) 

where the Lame constants are treated as invariant scalars in the transformation. The elastic relations in Eq. (|8(]) and 
( |Sl| ) are in the corotating frame S' . They can be obtained from the linear elastic relations in the inertial frame S, 
Eq. (|l9|), by using the transformations to the corotating frame in Eq. (36) and (^0|), and the coordinate transformation 
in Eq. (p5[). In doing the transformation of the elastic relations to the corotating frame, I am assuming that the tensor 
that enters in the inertial frame elastic relations in Eq. ( |l9| ) is the Eulerian strain tensor given in Eq.(20) and the 
quadratic terms have not been dropped. (As discussed earlier, dropping the quadratic terms in the displacement 
gradients is done in the corotating frame.) 

For small displacements in corotating frame S' , the gradients of the displacement vector in frame S' can be assumed 
small — for moderate angular velocity oj , so the quadratic terms in the gradient of the displacement can be neglected. 
Therefore, in the corotating frame S' , I take the relation between the radial component of the displacement vector, 
u, and the physical components of strain, e rr and e<p<p, to be (compare with Eq. (28) and (B9|)) 



du 

dr 1 

u 

- 



(84) 
(85) 



where the tilde on u indicates that the radial component of displacement field is taken in the corotating frame S' and 
the prime on r' indicates that the cylindrical radial coordinate is in the corotating frame S", see the transformation 
inEq. ©. 

Substituting Eq. ( p4) ) and ( pa ) into Eq. ( pp[ ) and ( pi] ) , leads to relations between the physical components of stress 
and the radial displacement field 



E 



du 
dr' 



u 



E 



u du 
v \ v o r 



(86) 
(87) 



Substituting these relations into Eq. (|7q ) leads to a differential equation for the displacement in the rotating frame S' 



,2 d u , du 



dr' 2 dr 



r — u 



E~ 



P u>y 3 



The general solution is 



U= E 



r 8 



(89) 
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Substitution of this solution into Eq. (|S6|)and ( |87| ) I obtain 

~° rr = C 1+ C 2; L-^p^r' 2 (90) 
*** = C,-C 2 ±- 2 - l -±^ P ^ (91) 

The stresses at r' = must remain finite, so I take Ci = 0. Applying the boundary condition on the long peripheral 
surface, Eq. ( |7^ ) leads to C\ = (3 + pw 2 6 2 /8 and the stresses 

cr- = ^p W 2 (fe 2 -r' 2 ) (92) 

= ip W 2 [(3 + ^& 2 -(l + 3^r' 2 ] (93) 
8 

The physical stress components in Eq. ( p2[ ) and (^) are in the corotating frame S'. However, due to the transformation 
between the corotating frame and the inertial frame in Eq. (|36|), and the coordinate transformation in Eq. (|35|), the 
corotating frame components in Eq. (|92] ) and (p3|) are equal to the inertial frame components of stress. Using the 
expressions in the rotating frame, such as Eq. (|75|) — (f77|), expressions for plane strain and other boundary conditions 
can be derived for rotating cylinders, disks and annular rings, see Ref. |], [|, ||, [| |). 



VI. SUMMARY 



The classic problem of stress in rotating disks or cylinders is important in applications to turbines, generators, and 
whenever large rotational speeds exist. The textbook problem of stress in perfectly elastic disks or cylinders is solved 
in standard texts |]J ||, ||, [|, ||, ||. The 'standard method' of solution begins with Eq. (|l|) and drops terms that are 
quadratic in strain gradient in the definition of the strain, see Eq. (|2p| ) . Equation (|l|) is valid only in an inertial frame 
of reference, since it is derived from Newton's second law of motion, which itself is only valid in an inertial reference 
frame. 

In this work, I have shown that dropping the terms quadratic in the displacement gradient (in Eq. (EOh) is incorrect 
in the inertial frame in which Eq. (|l|) is applied in the 'standard method' of solution jj], ||, |[ g, ^7 I provide an 
alternative formulation of the rotating elastic cylinder problem in a frame of reference that is corotating with the 
cylinder. In this corotating frame, I derive the dynamical equation for the stress (see Eq. ( |6l| ) or (|62|)) and I show that 
terms quadratic in the displacement gradient can be dropped because they are small (for moderate angular speed of 
rotation). This analysis in the corotating frame shows that the 'standard method' of solution jl], |3| |], g, g[ should 
be interpreted as being carried out in the corotating frame of reference of the cylinder. 

Furthermore, when stresses are computed in rotating disks or cylinders composed of materials that have more 
complex constitutive equations, such as elastic-plastic or viscoelastic behavior, one must carefully justify dropping 
the quadratic terms in displacement gradients. If dr opp ing these terms cannot be justified, then the problem can be 
analyzed in a rotating frame, using the derived Eq. (|6l]) or (^2|). Another practical application of the stress Eq. ( |62] ) 
in the rotating frame is to study elastic waves in bodies during rotation, where coriolis effects may play a role. 
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APPENDIX A: CONVENTIONS 



I specify tensor components on coordinate (non-holonomic) basis vectors using numerical indices, 1,2,3, such as 
a 12 . For physical components, which have the dimensions associated with that quantity, I use lettered indices, such 
as a r ^. In addition, I must distinguish between four coordinate systems: Cartesian and cylindrical coordinates in the 
inertial frame 5* and Cartesian and cylindrical in the corotating frame 5". I use z k — (x,y,z) and x k — (r,(j),z) for 
Cartesian and cylindrical coordinates in inertial frame S, respectively. In corotating frame <S", I use z' k = (x',y',z') 
and x' k = (r',cf>',z') for Cartesian and cylindrical coordinates, respectively. For distinguishing components in these 
four coordinate systems, I use an additional mark as follows: absence of mark and a bar, for Cartesian and cylindrical 
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TABLE I: Coordinates 

Inertial Frame S Rotating Frame S' 

Cartesian z k = z k z' k = z' k 

Cylindrical x k x' k 

TABLE II: Tensor Components 
Inertial Frame S Rotating Frame S' 

Cartesian a ab (z k ) = a ab , e lk , ui b a' ab (z' k ) = a' ab , e[ k , u' a » 
Cylindrical a ab (x k ), e lk , ui b a ab (x' k ), e lk , u) a h 



components in inertial frame 5*, respectively. For components in the corotating frame S' , I use a prime and a tilde, 
for Cartesian and cylindrical components, respectively. See Table I and II. 
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